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Abstract 

In a two species reaction diffusion system,we show that it is possible 
to generate a set of wavelength doubling bifuractions leading to spatially 
chaotic state. The wavelength doubling bifurcations are preceded by a sym- 
metry breaking transition which acts as a precursor. 
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Reaction Diffusion systems exhibit a rich variety of spatial patterns. The sim- 
plest such system consists of two species, one of which is autocatalytic and the 
other a fast diffusing antagonist. The interplay between the two gives rise to in- 
tersting Turing patterns [Q.A popular model is that of Gierer and Meinhardt 
|, g.It has the form 

A = DV^A + A^/B -A + a 

B = V^B + ^i{A^ - B^) 

Here A and B are the two interacting species. The three parameters in the prob- 
lem are D, the ratio of the diffusion coefficient of species A and B, a, the rate 
of production of A and fj, the reaction rate for the production of B.The ho- 
mogenous steady state A-^ = B = (1 + ct)^ is unstable to a patterned state 
along the boundary HcD = [(2/(1 -I- cr))^^^ — 1]^ and for fj, < for a given 
D. The homogenous steady state is unstable to a homogenous but temporal 
periodic state for /i < (1 — (t)/(1 + a).ln the region close to D ~ and /i ~ 1 
there is a chaotic state where the Hopf and Turing regions overlap |5[ |^. Steady 
but spatially chaotic states are however not known in the Gierer-Meinhardt 
model. On the other hand spatially chaotic states ^ are known to arise in the 
coupled map lattices ||, ^ |l^ which have recieved a lot of attention in the last 
decade. In this work, we have taken a two species model which is somewhat sim- 
pler than the Geirer- Meinhardt model above and explicitly shown the existance 
of a spatially chaotic state obtained by a process of wavelength doubling. It is 
interesting to note that the wavelength doubling is preceded by a symmetry 
breaking bifurcation similar to what occurs in a damped driven anharmonic 
harmonic oscillator [|ll|, |l2j. 

We consider a two species model where one of the species can be self satu- 
rating or self catalysing with a similar option for the other. The species interact 
via a linear cross coupling. Each species can help or hinder the growth of the 
other. The model for the populations A^i, A'2 of the two species is 

iVi = D\7^Ni + aiNi{l - iV|) + ^aA^i (1) 

N2 = V^A^a + 02^^2(1 - Nf) + h2Ni (2) 



2 



A spatially homogenous solution is one where both species die out i.e Ni = 
N2 = 0. The linear stability of the state in the absence of diffusion term shows 
that this will be a stable state of affairs if ai + 02 < and 0102 — 6162 > 0. We 
will consider oi, 02 , 61 and 62 to always satisfy thsc constraints. Consequently 
our focus will be on an instability brought about by the diffusion terms. The 
linear stability of Ni = N2 = against a perturbation SNi^2 = ^1,26*'^^ leads to 
the linearized equations 

di = {-Dk^ + ai)5i + biS2 (3) 

S2 = (-fc^ + 02)^2 + &2^1 (4) 
The growth rate p of the perturbation 61^2 can be found from 

Detf«i-^f-^ )=0 (5) 

\^ -62 a2-k^ - pj ^ ' 

so that 

2p = ai + 02 - (£> + l)fc^ ± 

+ 02 - (£> + l)fc2)2 + 4(aia2 - 6162 - k'^{Da2 + ai) + Dk^) (6) 

and instability can occur if 

aia2-hih2-k'^{Da2+ai) + Dk^ <Q (7) 

The extremum of the function in the left hand side of Eq(7) occurs at 

e = kl = \[a2+^-^) (8) 

and enforcing Eq(7) at this wavenumber, we need 

{a2D + aif > AD{aia2 - 6162) (9) 

This coiidition can be exactly satisfied and this implies that in the presence of 
diffusion ,the previously stable state can be destabilized in favour of a spatially 
periodic state. 

If we want to explore the stationary spatially periodic state ,then we need 
to write 

-/Vi,2 = Ai^2coskox (10) 

and find the coefficients Ai^2-This can be done by using the method of harmonic 

balance on Eq(l,2).Wc note immediately that terms of the form cosZk{)X will be 
generated. For the moment we ignore them and straightforward algebra leads to 

{ai-Dkl)Ar-^arA\ + brA2 = Q (11) 
(02 - kl)A2 - ^a2Al + 62^1 = (12) 
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We have verified that the cubic equation for that results from the above 
gives a unique answer i.e that there is only one positive root for a large range 
of values of the diffusion constant D. 

Returning now to the terms dropped in our harmonic balance we note that 
cos^kox produces a cosikox term which implies that the response will have a 
third harmonic as well, i.e the solution will be 

A''i^2 = Ai^2Cosk{)X + Bi 2Cos'6k(iX + . . . (13) 

It is easy to see that will now yield terms like Af 2Cos'^koxBi^2Cos3koX 
which will produce the fifth harmonic. Consequently in the expansion above all 
the odd harmonics will be present but not the even harmonics. This immediately 
means that the solution has a symmetry-namely if a; ^ a; + {2m + 1)^ then 
N{x + (2m + = —N{x). So long as a solution has this symmetry there 

cannot be a wavelength doubling bifurcation. 

We then ask the question if the above symmetry can be broken by the 
introduction of a cos2kox term. 

We try the stability of the solution of odd harmonics represented by Eq(13) 
by adding a perturbation SNi^2 = di^2e^*cos2kox. Straightforward algebra leads 
to the following qudratic equation for p, 

p^-{ai+a2-4.kl{D + l))p-4(aia2-bib2) + mkliDa2+ai)-16Dk^ = (14) 
where 

ai,2 = ai,2[l - ^{Al2 + BI2 + ^1,2^1,2)] (15) 

With ai + a2 < and q;iq;2 — &162 > it is clear that p can only be 
negative. In the presence of diffusion, instability will set in if 

Ak^{Da2 + ai) - 16Dk^ > aia2 - hb2 (16) 

We can tune the cross coupling coefficient or the self coupling to ensure that 
the above condition is satisfied. The broken-symmetry solution will look like 

-^1,2 = Ai^2coskox + Bi^2cos3kox + Ci^2cos2kox + . . . (17) 

Just like the coefficients Ai^2 and iJi,2 ,the coefficients Ci,2 can be found 
from harmonic balance to be given by 

(ai - Ak^,D)C, - ar(^ + ^ + + ^)ci = -61C2 (18) 

(«2 - OC2 - a2(^ + ^ + ^ + ^)C2 = -b2C, (19) 

The symmetry breaking solution of Eq(l,2) now allows wavelength doubling 
to take place. To see this we write 

A^i,2 = Ai^2coskox + Bi^2cos3kox + Ci^2cos2kox + 5i,2Cos^ (20) 
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Linearizing Eq(l,2) in 61^2 ,we find 

51 = -^61 + aiSi + biS2 - 3ai{Ai + Bi)Ci5i (21) 

52 = —^^2 + a2(52 + b25i - 3a2(^2 + -B2)C2^2 (22) 

By defining 

P1.2 = ai,2(l - 3(Ai,2 + -Bi,2)Ci,2) (23) 
and we can write Eqs(21) and Eqs(22) as 

5i = {Pi-^)5i+bi52 (24) 

4 = b25i + (/32 - ^)52 (25) 

The solution for 61^2 will be of the form e^'*. The condition for genration 
of cos^^ - the wavelength doubled state-can now be easily obtained in a form 
similar to Eq(16).0ur new state with anharmonic response can be written as 

-^1,2 = Ai^2coskox + Bi^2CX>s3kox + Ci^2cos2kox + Si^2Cos^^^ + ■ ■ . (26) 

where the amplitudes 81^2 can be found by using harmonic balance. To generate 
the next wavelength doubling one tries 

Ni,2 = Ai^2coskox + Bi,2Cos'Skox + Ci,2COs2kQX + Si,2Cos— — |-5i,2Cos— ^ (27) 

and linearizes Eq(l,2) in (5i^2-The resulting equations are 

ni.2 o 
-^'i = {ai - -j^)6i + biS2 - -aiAiSiSi (28) 

S2 = (a2 - -j^)S2 + 62<5i - -a2A2S2S2 (29) 

As we have done several times before,we can determine the condition for 
growth rate to be positive and thus the period quadrupled state is generated. The 
new state has the structure 

Ni^2 = Ai^2Coskox+Bi^2Cos'ikox+Ci^2Cos2kox+Si^2Cos—^+S'i 2(^08—^ (30) 

The sequence is now clear. We can add a perturbation ^i,2Cos^ and the term 
81^281 2 will trigger the instability that will lead to the generation of this longer 
wavelength state. A continous sequence of doubling can now take place leading 
to a state with no definite wavelengths. 

Thus we see that in a two-species model with one of the species autocatalytic 
and the other showing self saturation, we can have a spatially chaotic state at- 
tained through a succession of wavelength doubling bifurcations.lt is interesting 
that the wavelength doubling bifurcations has to be preceeded by a symmetry 
breaking bifurcation. 
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